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Abstract. The intensity of the soft bremsstrahlung depends on the colour 
charges of the hard- scattered partons. This intensity is larger in dijet production 
with resolved photons by the factor Ca/Cf = 9/4 in comparison with direct 
photoproduction for particular configurations. It is investigated how these effects 
are reflected in Monte Carlo models. 



The idea that the shape of inclusive hadron spectra can be derived from the 
perturbatively calculated parton distributions (Local Parton Hadron Duality Q) has 
been quite successful in general. There are some specific predictions for the very soft 
particles. Because of their long wavelength the soft gluons "probe" only the colour 
factors of the primary hard partons (the "colour antennae patterns" (jjj) and then 
their production properties can be derived directly from the lowest order perturbative 
diagrams. Recently the case of soft particles has been investigated in more detail 
||. In the simplest case of e + e _ — > qq the density of the soft gluons is predicted to 
be practically energy independent. Remarkably, it was found that the soft hadrons 
indeed follow this expectation suggesting the relevance of the duality picture even for 
low momentum particles (p <~ 300 MeV). 

This success could be accidental and may reflect purely hadronic properties. It 
is therefore important to test more specifically the perturbative origin of the soft 
phenomena. A direct consequence of the dominance of the lowest order diagrams is 
the proportionality of the soft particle density to the primary parton colour factor 
which is Cf = § and Ca = 3 for a qq and gg antenna respectively. Then low 
momentum particles are produced with the relative intensity |. It is important to 
recall that the total event multiplicities approach this ratio rather slowly with energy. 
As a gg system is difficult to prepare one has to find an equivalent parton antenna 
with the appropriate effective colour factors. Such a realization is provided by, for 
example, dijet photoproduction. 



1. Tests with Photoproduction of Dijets 

One can distinguish in the leading order QCD approach the direct and the resolved 
processes. In the first case the photon participates directly in the hard scattering 
subprocess by photon-gluon-fusion or QCD-Compton scattering and transfers a large 
fraction (x 7 ~ 1) of its primary energy to the secondary jets. In the second case, the 
hard scattering subprocess involves the partons (q, q and g) in the photon and in the 
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proton and the energy fraction x 1 < 1. At HERA the dominant direct and resolved 
processes correspond to quark and gluon exchange respectively and the expected 
distributions in the dijet cms scattering angle S have been clearly observed at HERA 

!■ 

Particularly simple results arise for the soft radiation perpendicular to the 
scattering plane In this case complications with the cut-off in the transverse 
momentum fc_L > Qa disappear and all formulae depend only on the angles between 
the hard partons. In the simplest case of a qq dipole the density of gluons with 
momentum p perpendicular to qq is 

The soft radiation in more complicated hard processes involving gluons, relevant in 
the present discussion, can be obtained from appropriate superpositions of elementary 
dipoles. 

We consider the ratio R\ of the perpendicular radiation in the process i and in 
a standard qq process with O q q = tt. With this normalization one finds for the direct 
processes a) 717 — > qq and b) 717 — » qg 



R c \ =1, R\ V 



1 



3 — cos 6 S — ^Tj (1 + cos ©s 



(2) 



In both cases R l ± — > 1 for scattering angle <d s — > 0, i.e. in this limit both processes 
behave like qq antennae. In the QCD-Compton process b) the qq antenna changes 
into a gg type antenna for O s — ► n. 

In case of the dominant resolved processes the results in Eqs. (58-61) were derived 
[0 for small scattering angles. There is only a weak dependence on this angle and for 
O s — > 0, where the gluon exchange dominates, they approach 

R\ = Ca/C f = 2.25. (3) 

Using the known full expressions for the soft gluon radiation patterns B we found the 
approximate results for the leading processes gg — > gg and qg — > qg to be correct within 
~ 10% for O s < 7r/2; within this accuracy one can safely neglect the contribution 
from the process gg — » qq (see Eqs. (A7-A9) in [4c]). In the table we show the 
results for the ratios R± in the case when the final partons are not identified, i.e. 
after symmetrization as in Eq. (53) of 0. One can see that the result (|j) remains 
approximately valid at arbitrary scattering angles within ~ 20%. 



cos Q s 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0.0 

73 -> qq 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 1.0 

79 -> qg 1.0 1.16 1.29 1.39 1.46 1.52 1.56 1.59 1.61 1.62 1.62 

gg -> gg 2.25 2.31 2.36 2.41 2.46 2.51 2.55 2.58 2.61 2.62 2.62 

qg -► qg 2.25 2.23 2.19 2.15 2.10 2.04 1.98 1.92 1.87 1.84 1.83 



Table 1. Angular dependence of symmetrized cross section ratios R± for quark 
and leading gluon exchange processes 
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2. Monte Carlo Studies 

We have made a 'Monte Carlo measurement' of this ratio using the HERWIG § 
program. HERWIG contains coherent QCD radiation effects and a cluster 
hadronization model. However, since HERWIG produces a complete simulation of 
events, realistic jet algorithms and detector acceptance cuts can be applied. Thus, 
if HERWIG is consistent with the analytical results, this is an indication that the 
perturbative calculations are relevant for soft particles as well, that the results are 
insensitive to detector acceptance effects, and that the measurement is likely to be 
feasible. 

HERWIG events with a hard scatter of px > 6 GeV were generated for HERA 
beam conditions in the region 0.2 < y < 0.85, Q 2 < 1 GeV 2 , where y is the usual 
inelasticity variable and Q 2 is the virtuality of the photon exchanged between the 
proton and the positron. The Ut jet finder was run in inclusive mode [|| on the 
hadronic final state, and two jets with E T ct > 6, at least one of which must have 
E T et > 7 GeV, were demanded, in the pseudorapidity region \rj\ < 2. We have also 
imposed the cut \rji+r)2\/2 < 1, constraining the boost of the dijet system so to achieve 
a more uniform acceptance in jet scattering angle, as described in ||. Together, these 
cuts correspond to a region in which measurements have previously be made at HERA. 
Only 4.5 pb _1 of simulated data was used. 

The particle pt spectrum down to 50 MeV for particles within a cone of one unit 
in i] — 4> from the vector perpendicular to the dijet system in the dijet centre of mass 
system is then obtained. 

Based upon the variable x® BS Q , the fraction of the photon's momentum which 
enters into the dijet system, we then divide the events into two samples - i° BS > 0.75 
('resolved') and x° BS < 0.75 ('direct'). The ratio of the pr spectra for resolved/direct 
events is calculated, and is shown in the upper two figures for two regions of scattering 
angle, large (left) and small (right) - solid points. The small scattering angle region 
is defined by \Arj\ > 2 and the large by \Arj\ < 2. In terms of cos6 s , these ranges 
correspond approximately to < cos 8., < 0.76 and 0.76 < cos0 2 < 0.96. The 
plots also show the predictions above for these kinematic ranges. The line shows the 
prediction taking the subprocess mixture as given by HERWIG. The band shows the 
uncertainty in the prediction if absolutely no knowledge of the partonic subprocess 
type, but a perfect separation of resolved and direct type diagrams is assumed. In all 
cases, processes not included in the table are neglected. 

The HERWIG results approach the analytic predictions at low px- Also shown 
(clear circles) is the result from HERWIG when multiparton interactions are allowed |jj 
with ap™" for the hard scatter of 2.5 GeV. Multiparton interactions are often appealed 
to as a means of improving agreement between data and MC in the HERA region. 
They raise the ratio slightly, moving it closer to the prediction for high angle scattering, 
and further away for small angle scattering. 

Because the resolved cross section peaks strongly at low scattering angles while 
the direct matrix element does not, the high x® BS sample for small scattering angles 
in fact consists of 30% resolved-type diagrams. This is why the 'prediction' lies outside 
the band, since in constructing the band it was assumed that the separation between 
resolved and direct LO diagrams was perfect. The extra contamination from resolved 
acts against the enhancement in the ratio coming from the angular dependence of the 
radiation. Clearly a better way of distinguishing resolved and direct type diagrams is 
required - either a new variable, or a harder cut on i° BS . We have tried cutting at 
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Figure 1. Ratio of the average track transverse momentum in high s to low 
3.0BS everi i s 'Resolved/ 'Direct') for large jet scattering angles (left) and for 
small jet scattering angles (right). 



= 0.9, at the same moving the separation between low and high angular regions 
to | A?/ 1 = 1.5(cos8 s = 0.64). The results are shown in the lower plots. The prediction 
for low angle scattering is raised slightly and the agreement with HERWIG is rather 
good. So it looks like for these distributions the HERWIG results reproduce duality 
picture quite closely. Further improvements in these tests are undoubtedly possible. 

The results presented reiterate that it is of great interest to measure these 
distributions in HERA data. The direct comparison of the low momentum particle 
production in direct and resolved processes should yield the ratio of colour factors 
Ca/Cf- Its observation is a crucial test of the significance (or otherwise) of 
multiparton interaction models and, more fundamentally, of the perturbative picture 
for soft particle production. 
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